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1. Introduction 

The idea tliat our Universe is a defect embedded in a higher- dimensional bulk space [T] 
provides a context in which the cosmological evolution can display novel features. Of 
particular interest is the possibility of exchange of energy between the defect and the 
bulk space. The first obstacle in realizing this idea is the requirement to localize low- 
energy gravity on the defect. This is necessary in order to reproduce the correct form 
of Newton's law, and conventional cosmological evolution with the Hubble parameter 
depending linearly on the energy density of the matter localized on the defect. The 
Randall- Sundrum model |2I provides an example in which this requirement is satisfied. 
The defect is a four-dimensional hypersurface (a 3-brane) in a five-dimensional bulk with 
negative cosmological constant (AdS space). The geometry is non-trivial (warped) along 
the fourth spatial dimension, so that an effective localization of low-energy gravity takes 
place near the brane. The matter is assumed to be concentrated only on the brane. 
For low matter densities the cosmological evolution on the brane becomes typical of a 
Friedmann- Robertson- Walker Universe jHI E] . 

It is possible to discuss the brane evolution either in a coordinate system (system A) 
in which the brane is located at a fixed value of the fourth spatial coordinate and the bulk 
is time-dependent 13] , or in a different coordinate system (system B) in which the bulk 
is static and the brane is moving 0. In the latter case, the bulk metric is five-dimensional 
AdS-Schwarzschild jH]. The two points of view are equivalent [7j. Several aspects of the 
cosmological evolution have been the subject of a large number of studies. (For recent 
reviews, with extensive lists of references, see ref. j^.) 

We are interested in the possibility of energy exchange between the brane and the 
bulk, and the modifications it would induce on the cosmological evolution on the brane. 
For such a study the bulk must be assumed to contain, apart from the cosmological 
constant, some matter component. The basic equations have been studied in the 
coordinate system A. However, in this system it is difficult to obtain an exact solution for 
the bulk metric. In the presence of bulk matter, the cosmological evolution on the brane 
is not autonomous, and the explicit knowledge of the bulk energy-momentum tensor is 
necessary. As a result, certain assumptions often need to be made for its form P 
Previous studies have considered a bulk field ^2j, or radiation emitted by the brane into 
the bulk ^3 Also an attempt has been made to incorporate thermal effects in the 
bulk US]. 

In this study we start by deriving an explicit solution for the bulk in the coordinate 
system B. We allow for a negative cosmological constant and additional matter with an 
arbitrary equation of state. In the following section we obtain a particular class of static 
solutions, that correspond to generalizations of the AdS-Schwarzschild metric. These 
solutions can be interpreted as AdS-stars, as their form resembles the solution for the 
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interior of stellar configurations. In section 3 wc introduce a brane as the boundary of the 
five-dimensional space. We employ the coordinate system A, in which the brane is located 
at a fixed value of the fourth spatial coordinate. We also establish the connection between 
the two coordinate systems. This permits us to derive the exact equations governing the 
cosmological evolution on the brane, making use of the exphcit solutions of section 2. 
In section 4 we study the cosmological evolution on the brane. The solutions describe 
energy outflow from the brane and the presence of "mirage" contributions to the Hubble 
expansion. In the final section we summarize our results. 



2. The bulk 



We consider an action of the form 

j (fx ^ {M^R + A + + J d^x^/^ {-V + Z:™"*) , (2.1) 

where R is the curvature scalar of the five-dimcnsional metric qab, A, B = 0, 1, 2, 3, 4, —A 
the bulk cosmological constant (A > 0), jab, ^ — 0,1,2,3, the induced metric on the 
3-brane, and V the brane tension. 

We assume that the metric in the bulk can be put in the static form 

ds^ = -n\r)dt^ + r^dnl + b''{r)dr\ (2.2) 

where dQ^ is the metric in a maximally symmetric three-dimensional space. We use 
k = —1, 0, 1 to parametrize the spatial curvature. The non-zero components of the five- 
dimensional Einstein tensor are 
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where T\ denotes the total energy-momentum tensor. Assuming a negative cosmological 
constant —A and a static perfect fluid in the bulk, we write the bulk energy-momentum 
tensor as 



T\ = diag(A - p, A p, A p, A p, A p). 
Combining the 00 and 44 Einstein equations we obtain 



(6n)' 
bn 
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/^From the conservation of the energy-momentum tensor we have 

- = (2.9) 

n P + P 

In order to solve the Einstein equations we proceed in complete analogy to the 
standard study of stellar structure JB]- We rewrite the 00 Einstein equation as 

^) =2^- + ^-'(A-p)- (2.10) 
This has the solution 

where Ai{r) satisfies 

^ = 2vrVp. (2.12) 
dr 

For p = and constant Ai we obtain the known AdS-Schwarzschild solution jB] . Through 
the use of eqs. ()2.9p . ()2.11|) . the 44 Einstein equation becomes 



p' 11 



-1 



p + p r r 

The solution of eqs. ()2.9|) . ()2.1H) - ()2.13|) automatically satisfies eqs. fl2.4|) . ()2.8|) that were 
not used in the derivation. 

For a complete solution the equation of state p = p{p) must be specified. Then, 
if initial conditions are given for p and Ai at some r, the system of equations ()2.12|) . 
fl2.1H|l can be integrated. The resulting solution determines b{r) through substitution in 
eq. ()2.11|) . Finally, the integration of eq. ()2.8|) or ()2.9|) determines n(r). 

The existence of a solution depends crucially on our assumption for the form of the 
metric. In the absence of matter in the bulk, the most general solution is the AdS- 
Schwarzschild that can be put in the form of the ansatz ()2.2|) . This solutions results 
from our equations under the assumption p = p = 0: Eq. ()2.8j) gives b = 1/n, while the 
integration of eq. ()2.10p gives eq. 1)2.111) with constant Ai. 

In the presence of matter in the bulk, the ansatz of eq. ()2.2j) is not the most general. 
However, as we shall show in the following, our assumption leads to a physically motivated 
solution, that provides an interesting example of cosmological evolution induced on the 
brane. Physically, the distribution of matter in the bulk corresponds to that of a stellar 
object in AdS space. Most of the matter is concentrated near the center of the star at 
r = 0, but a non-zero density exists for large r as well. In this sense, the brane moves 
within a background that can be characterized as an AdS-star configuration. If all the 
matter falls into a black- hole singularity at r = 0, our solution reduces to the standard 
AdS-Schwarzschild metric in five dimensions. 



.(2.13) 
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In a pure AdS background, described by eq. with n^(r) = l/6^(r) = 

k + Ar^/(12M^), the concentration of matter is expected to be enhanced near r = 0. 
We are interested in solutions with finite energy density throughout the bulk. For this 
reason we consider the boundary condition p(r = 0) = po, p{r = 0) = Pq. The solution 
of eq. ()2.13|) is now uniquely specified. For the solution of eq. ()2.12p we assume that 
Ai{r = 0) =0. In this way Ai{r) can be roughly interpreted as the integrated matter 
between r = and some arbitrary value of r. In order to prevent the last term in brackets 
in the r.h.s. of eq. ()2.1H|1 from vanishing at a certain r, we are forced to take k = 1. 
This implies that the brane Universe must have positive spatial curvature. The functions 
l/6^(r) and n'^{r) never vanish. The resulting solution is finite everywhere in the bulk 
and has no horizons. For the cases k = 0,-1 the function l/6^(r) vanishes at a certain 
r, which forces n^(r) to vanish at the same point, as can be seen from eq. ()2.8|) . As a 
result, our ansatz assumes a static distribution of matter in the presence of horizons in 
the metric. This requirement cannot be satisfied, as the divergence of p implied by eq. 
()2.13|1 indicates. It is probable that an alternative non-static ansatz for the metric (such 
as the generalization of the Vaidya metric [Hi) must be considered. 

The solution of eqs. ()2.12j) . ()2.13p with k = 1 has strong similarities with the well 
known solutions for the interior of stellar configurations. We cannot provide an analytical 
solution for all r, but we check our conclusions through the numerical integration of eqs. 
()2.12|1 . ()2.13|1 . For a generic polytropic equation of state p = wp^ and the boundary 
conditions we chose, the solution is regular for small r with p = pq, A4{r) ~ Por'^-t 
The large-r behaviour can be obtained analytically through the assumption that the 
cosmological constant term ~ Ar^ dominates over the curvature and matter terms. In 
this way we obtain 

1 

/ 1-7 1 \ 7-1 

p(r) = [ nr for77^1 (2.14) 

. w-\-l , , 

p(r) = nr » for 7 = 1, (2.15) 

with K an integration constant. Regular behaviour near r = is always connected with 
this type of asymptotic behaviour, as we have verified numerically. 

Let us consider first the large-r behaviour for 7 7^ 1. For 7 < 1 (we always assume 
7 > 0), we obtain p ~ r^^/"^, p ~ r~^. The density of bulk matter is non-zero for all r. 
We obtain 

llTl^ 1 47-1 

M{r) = — K— r— + M^, (2.16) 

47 — 1 

with another integration constant. For 7 > 1 there is a finite value of for which 
p = p = 0. At this point eq. ()2.8|1 implies that b{rs) = l/n{rs). We can assume that 

I For 7 = 1 there is another solution that diverges near r = according to p ^ r^^, A4 ^ r^. This is 
similar to the standard neutron star solution jTB| . 
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the matter density remains zero for r > r^. This permits the matching of the solution for 
r < rg with the standard AdS-Schwarzschild metric for r > r^. Simply, in the latter we 
must take A4 = Ai{rs). This is analogous to the matching of the interior solution with 
the Schwarzschild metric at the surface of a star. 

In cosmological applications the linear equation of state with 7 = 1 is most commonly 
employed. For large r and w > the matter density never becomes exactly zero, even 
though it falls off as a power of r. We obtain 

, , , , 2w7r^ 31^-1 , ^ , „x 

M{r) = -nr— + M^, 2.17 

6W — 1 

with an integration constant. For the particular value w = 1/3 the asymptotic 
behaviour is 

Mir) =27r\ln(^^^ +M^, (2.18) 
with Tfe another integration constant. 

3. The brane 

For the discussion of the cosmological evolution on the brane we follow ref. jSj. First we 
consider a coordinate system (using Gaussian normal coordinates) in which the metric 
takes the form 

ds"^ = 'jabdx'^dx'' + drj^ = -m^(f, ri)df'^ + a^(f, r])dill + drj^ . (3.1) 

The brane is located at = 0. In order to keep our discussion simple we identify the 
half-space 77 > with the half-space < 0, in complete analogy to ref. [2]. We define the 
proper time on the brane dr = m(f, rj = 0) df. We can rewrite the above metric as 

ds'^ = -m^{r, ri)dT^ + a^{T, ri)dnl + drj^, (3.2) 

with m(r, rj) = m(f, ri)/m{f, rj = 0), a(r, rj) = a(f, rj), so that m(r, 77 = 0) = 1. 
Through an appropriate coordinate transformation 

t = tir,r]), r = r{T,r]) (3.3) 

the metric (j2.2j) can be written in the form of eq. (j3.2j) . It is obvious that r = a(r, 77). 
Therefore, we obtain 

9t , dt ^ 

dt = —dr + —dr] (3.4) 
OT or] 

dv dv d(x 
dr = T^dr + —dr] = T^dr + —dr]- (3-5) 
OT or} or or] 
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We define -R(t) = a(r, 77 = 0). In the system of coordinates (t,r) of eq. ()2.2|1 the 
brane is moving, as it is located at r = R{t). At the location of the brane, we find 



dt 



dr n{R) 
dt h{R) 
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(3.9) 

(3.10) 

(3.11) 
(3.12) 



where the dot denotes a derivative with respect to proper time. The negative sign in the 
r.h.s. of eq. ()3.12j) is a consequence of our assumption that r decreases on both sides of 
the brane 

For completence, we also give the components of the brane velocity u in the coordinate 
system (t, r): 



u 



u 



n{R) 
R. 



b\R)R^ 



nl/2 



The unit vector n normal to the boundary defined by the brane motion has 
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n{R) 
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(3.13) 
(3.14) 

(3.15) 
(3.16) 



It is now straightforward to express the bulk energy-momentum tensor at the location 
of the brane in the coordinate system (r, r^). We find 

T° = A - p{R) - \p{R) + p{R)\ b'^{R)R^ (3.17) 

T\ = T\ = T% = A + p{R) (3.18) 
T\ =A + p{R)+\p{R)+p{R)\b\R)R' (3.19) 



1/2 



b{R)R b\R)R' + l p{R)+p{R) 



(3.20) 
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We have indicated explicitly the dependence of p, p and 6 on i? in order to emphasize 
that they are the solutions of eqs. ()2.11|) - ()2.13|) . with R as their argument. 

The equations governing the cosmological evolution on the brane can be obtained in 
various ways. The most straightforward is to write the Einstein equations in the system 
(r, rj), assuming that the brane energy momentum tensor has the form 

Tl = 6{r]) diag{V - p,V + p,V + p,V + p,0), (3.21) 

with V the brane tension and p, p the energy density and pressure of a perfect fluid 
localized on the brane. The details of the calculation can be found in refs. [Hl^JI] and 
will not be repeated here. The evolution on the brane is governed by the equations 

P + 3|(p + p)= -2n (3.22) 

X + ^Hx =^(/^ + ^)^4 (3.24) 
+ 4i^0 = - {T\ - A) , (3.25) 

where A = {V'^/12M'^ — A)/12M^ is the effective cosmological constant, and T\, T\ are 
given by eqs. ()3.19|) . ()3.20|) . For the particular form of the energy- momentum tensor that 
we are considering, we find that 

satisfies eqs. (jSI21|), (jH^ . if M{R) is the solution of eq. (|TT^ with r replaced by R. In 
order to show this, we make use of eq. (|3.23|) in the form 

-H'- ^ (?? + 7Vn\ ^ + ^_M{R) 
rT--^ -TuaT^^P + 2^^) - ^ + 6vr2M3 R^ +^ 

^ iP + yf-U7^.- (3-27) 



144M6 ' b'^{R)R'^' 

We can now express eq. ()3.20|) as 

n = Y^s^^R) R' H {p + V) [piR) + piR)] . (3.28) 

The rate of energy trasfer is fixed by the static solution for the bulk that we derived in 
section 2. This solution is consistent only with energy loss by the brane. In order to 
understand better this point it is convenient to employ the coordinate system B as in 
section 2. In this system the brane forms the moving boundary of a static bulk. As the 
boundary shifts towards larger values of the coordinate r, energy must be transferred from 
the brane to the bulk. Our solution corresponds to the situation that matter leaves the 
boundary fast enough so as to fill the additional volume before flows develop from the rest 
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of the bulk space. Of course there are constraints on the reahzation of such a scenario. If 
the brane matter cannot escape fast enough the bulk cannot remain static. The detailed 
discussion of this issue must deal with complicated issues, such as the exact profile of the 
brane beyond the approximation by a (5-function, the mechanism that localizes particles 
on the brane, the dynamics near the surface (pressure gradients, local flows). In this work 
we simply assume that matter can escape from the brane fast enough for our solution to 
be realized. 

An alternative way to derive eq. ()3.27j] is by considering the junction conditions for 
the boundary defined by the brane motion. In the coordinates (r, rj) of eq. ()3.2|) the 
extrinsic curvature of the boundary is 

Kab = ^d,-fab. (3.29) 
The junction conditions are 

AK^b = Kb - Kb = (^^5 - Ir^lab^ , (3.30) 

where it denote the two sides of the brane, and the energy-momentum tensor is given by 
eq. (|3.21|) neglecting the (5-function. For the spatial components we find 

Then eq. fl3.30|) immediately reproduces eq. ()3.27|) . 



4. The cosmological evolution 

The cosmological evolution on the brane is described by eqs. ()3.22p . ()3.27p . with given 
by eq. or (nOHjl . and p{R), p{R), M{R), b{R) determined by the solution of eqs. 

fl2.11|) - ()2.13|) that we discussed in section 2. 

The evolution has some characteristic properties: The energy on the brane is diluted 
by the expansion, but is also reduced through the outflow from both sides of the brane. 
This implies that the energy density of the various components of brane matter (dust, 
radiation etc) does not scale as function of R with the "naive" power determined 
by the equation of state p = wp. The outflow makes the reduction faster. 

The effective Friedmann equation ()3.27|) for the expansion on the brane has several 
terms: The effective cosmological constant A leads to exponential expansion. In order to 
get conventional cosmology we have to set it to zero through the appropriate fine-tuning 
of the bulk cosmogical constant and the brane tension: V"^ = 12M^A. The curvature 
term has the standard form. In the context of the solution of section 2 for the bulk we 
have k = 1. The brane matter induces a term ~ p^, characteristic of the unconventional 
early cosmology of the Randall- Sundrum model 13]. At low energy densities p <^ V , 
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the standard contribution ~ p to the Hubble parameter dominates. The effective Planck 
constant is M^j = 12M^/V. 

The novel feature of eq. ()3.27p is the term ~ J\4{R)/R'^. It is a generalization of the 
"mirage" , or "Weyl" , or "dark" radiation term [SI CUl CSl UH US] encountered in brane 
cosmology, which is related to the mass of the black hole in the bulk. We recover the 
1/i?^ dependence, modulated by the explicit i?-dependence of the mass term JH. Similar 
modifications have been found in several variations of the basic Randall-Sundrum scenario 

(nnHiin]. 

For a bulk described by the solutions of section 2 the cosmological expansion on the 
brane at early times is dominated by the brane matter density. The reason is that the 
bulk effects are incorporated in the term ~ Ai (R) / R'^ that becomes constant in the limit 
i? — >• 0. The energy outflow T°4 = R^Hp{pQ + po) /(12M^) also becomes subleading in 
this limit, and the energy density on the brane is reduced mainly through the expansion: 
p ~ R~^^^^^\ The scale factor obeys R ~ ^i/(3+3«i>)^ similarly to the high-energy regime 
of the standard Randall-Sundrum model. 

We are mainly interested in the large- i? behaviour of the cosmological solution. In 
this limit, the form of p{R) has been determined analytically at the end of section 2, and 
is given by eqs. ()2.14j) . ()2.15|1 with r replaced by R. For 7 > 1 the energy density of 
bulk matter vanishes for a certain R and remains zero subsequently. The cosmological 
evolution on the brane has a complicated behaviour up to a certain time. Later on, 
it becomes conventional, with the "mirage" radiation term the only indication of the 
presence of the bulk. 

For 7 < 1 and large R we have p{R) ~ R~^^"', p{R) ~ -R~^-§ We also obtain 
1/6^ = Ai?2/(12M3) and T\ = HV{p + p)/A. In this way we find that for large R the 
energy density on the brane obeys d{p R^'^^^'^'') / dR ~ —R^^^"^. For w > this implies 
that the energy density becomes zero at some finite Rj. Our solution does not have a 
physical meaning beyond this point, as it describes matter with negative energy density. 
For Rf > R ^ 1 the form oiAi{R), as given by eq. ()2.16|) . implies that there is a "mirage" 
contribution to the Hubble parameter that scales ~ R~^^"', as well as a "mirage" radiation 
term. The former contribution dominates the radiation term for 7 > 1/4 and sufficiently 
large R. The influence of the bulk matter on the brane evolution is equivalent to that of 
a "mirage" matter component with a linear equation of state and Weff = — 1 + 1/(87). 
For 1/3<7<1, we have We// < 0. For 1/2 < 7 < 1 the "mirage" term dominates the 
curvature term for large R. 

The most common form of the equation of state in cosmological applications has 

7 = 1. This choice leads, through use of eq. ^J^, to d{p R^^^-^^^) / dR /?i+3u.-i/«;^ Yoi 

w < 1/ (2 + 3w) the energy density p never vanishes, but falls off asymptotically according 

§ We emphasize that, for 7^1, the scale factor R must be sufhciently large for the term ^ i?(i-7)/7 ^ 
dominate over 1/w in eq. (|2.14|) . 
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to the standard relation p ~ 

^-3{i+u,)_ pQ^y^ > 1/(2 + 3w) the energy density p becomes 
zero at a finite Rf, beyond which our solution does not have a physical meaning. Eqs. 
fl2.17|) . ()3.27|) imply that the infiuence of the bulk is equivalent to a "mirage" radiation 
term and another "mirage" component with Weff = —2/3 + l/{3w). For w > 1/2 we 
have Weff < 0. For a "stiff" equation of state with w = 1, the "mirage" term acts as an 
additional curvature term for large R. 

The most attractive from the above scenaria assumes linear equations of state for 
both the bulk and brane matter with w < 1/(2 + 3w). For large R the energy outfiow 
becomes subleading, and the brane energy density falls off according to the standard 
relation ~ p^^ w < 1/3 the main "mirage" effect is a radiation term. For 

w = 1/3 this term gets logaritmic corrections ~ ln{R)/R'^. For 1/3 < w < 1/(2 + 3w) 
the leading "mirage" effect is a component with an effective equation of state with 
Weff = —2/3 + l/{3w). In all the scenaria the "mirage" effects are negligible at early 
times, when the brane energy density dominates the evolution. However, at later times 
this energy is depleted through the expansion and the outflow and the "mirage" terms 
(as well as the curvature term) dominate the Hubble expansion. 

As a flnal remark we point out that in the limit of vanishing matter energy density 
both in the bulk and the brane our solution reduces to the Randall-Sundrum solution. 
As a result we expect localization of the massless Kaluza-Klein mode of the graviton 
around the brane. This is consistent with the emergence of a cosmological evolution 
typical of a four-dimensional Universe for brane energy densities smaller than the brane 
tension and bulk energy densities smaller than the absolute value of the cosmological 
constant. According to our discussion at the end of section 2, the cosmological constant 
term always dominates the bulk metric for large values of the scale factor R. Similarly, 
the brane tension dominates the brane energy density in the same limit. 

5. Conclusions 

The general picture that emerges from our discussion is that the presence of matter in the 
bulk affects the brane evolution in unexpected ways. The "mirage" radiation term, that 
has appeared in previous studies, is only the simplest manifestation of the existence of 
the bulk. It comes about because of the presence of a black hole at r = in an AdS bulk, 
in the static coordinate system of eq. ()2.2j) . In this study we constructed more general 
bulk solutions that permit a non-trivial distribution of bulk matter. Because of the bulk 
geometry most of this matter is concentrated in the region near r = in the coordinates 
of eq. ()2.2j) . or for 77 — oo in the coordinates of eq. ()3.2|) . However, it is possible for the 
brane to be located in a region of non-zero density of bulk matter. As a consequence, the 
cosmological evolution on the brane is modifled. 

The first important effect is the flow of energy away or towards the brane. The 
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solutions we constructed are consistent only with energy outflow, accompanied with a 
faster than usual dilution of brane matter. However, it is possible that more sophisticated 
bulk solutions could permit the influx of energy, accompanied with exciting novel effects 
[inilll]. The second consequence is the modification of the standard Hubble expansion 
through "mirage" contributions to the brane energy density. The effective equations of 
state of such components depend on the equation of state of the bulk matter in a non- 
linear way. For example, bulk matter with a linear equation of state and w > may lead 
to a "mirage" component on the brane with Weff < 0. A general property of the solutions 
we derived is that the "mirage" effects are negligible at early times when the brane energy 
density dominates, but become significant at late times. 

The bulk solutions we derived are only a small sub-class of the more general solutions 
that could predict interesting low energy cosmological evolution for a brane Universe. 
The novel features introduced by the presence of matter in the AdS bulk is the outflow 
or influx of energy from the brane, and the appearance of "mirage" components affecting 
the cosmological evolution on the brane. Such features could be useful for modifying 
the evolution so as to include, for example, an accelerating implied by recent 

astrophysical observations The task of deriving new bulk solutions is difficult, but 

the implications appear to be very interesting. 
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